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Abstract 

We study the volume dependence of the left-current correlator with non-degenerate quark masses to next- 
to-leading order in the chiral expansion. We consider three possible regimes: all quark masses are in the 
e-regime, all are in the p-regime and a mixed-regime where the lighest quark masses satisfy m^T,V < 1 
while the heavier nisTiV 3> 1. These results can be used to match lattice QCD and the Chiral Effective 
Theory in a large but finite box in which the Compton wavelength of the lightest pions is of the order of 
the box size. We consider both the full and partially-quenched results. 
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1. Introduction 



The first principles determination of the low-energy couplings of the Chiral Effective Theory, that de- 
scribes the meson interactions at low momenta, is one of the milestones of lattice QCD. This matching 
can only be carried out reliably close to the chiral limit, and this is often a limitation for lattice simula- 
tions, because the computational cost of lattice simulations increases very significantly with decreasing 
quark masses and increasing volumes. 

Even though there has been important algorithmic progress in recent years, it seems quite difficult to 
reach the range of the u and d quark masses, at least within the p-regime, i.e. satisfying the condition 
MttL » 1. 

The e-regimc [1, 2] has been advocated [3] as an alternative to perform the matching, that a priori could 
be more economic in the sense that the quark masses can be taken to zero without increasing the box-size 
proportionally, since in this regime M^L < 1. Finize-size effects are large in this situation, but they are 
calculable within the Chiral Effective Theory in terms of the infinite- volume low-energy constants [4]. 
Even though the truly chiral regime requires that the volume is scaled to infinity eventually, not just the 
quark mass to zero, the scaling with the volume at zero quark mass is more universal in the sense that 
it involves less low-energy couplings, since most of the operators that appear at higher orders in Chiral 
Perturbation Theory (ChPT) include explicit powers of the quark mass. 

In particular it can be shown that in the e-regime only the leading order couplings F and S appear in 
two-point functions^ at next-to-leading order (NLO) of ChPT. 

A number of correlation functions have been computed to NLO in the e-regime of ChPT. The two-point 
functions of scalar, pseudoscalar, vector and axial-vector correlators were presented in [4]. In [5, 6] the 
same correlators were computed in quenched ChPT and also in the full theory, but in a fixed topological 
sector [12]. Three-point functions relevant for determining the weak low-energy couplings were first 
presented in [8] both for the full and quenched theories. The e-regimes has also been recently applied 
to the study of baryon properties [9]. These results have been used in various simulations to extract 
low-energy couplings mostly in the quenched approximation [3], but more recently also in unquenched 
simulations [10]. 

In many of these computations quark masses were taken degenerate. The purpose of this paper is 
to present the results for the left-current correlator for non-degenerate quark masses. The methods 
developed here can be readily applied to other correlation functions. 

More concretely, we will consider the situation with Ng heavier quarks with common mass and iVy 
light quarks with common mass m^. At this point one could imagine having three different situations: 

• m^isTiV ^ 1 : all quarks are in the p-regime. 

• m^igTiV <C 1: all quarks are in the e-regime. 

• m^TiV < 1, rngTiV ^ 1: some quarks are in the e and some in the p regime. It appears natural to 
identify m„ with the light quarks u, d and with the s quark mass in realistic simulations. 

^Up to contact terms. 
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The explicit and Ng dependences will be shown in such a way that the partial quenching of the 
w or s quarks can be easily done via the replica method [11]. Considering the partial quenching of the 
lighter quarks is interesting in the context of mixed-actions [19], where the valence and sea quarks are 
treated in different regularizations, for example with overlap valence quarks and Wilson sea quarks. 

The structure of the paper is as follows. In section 2, we present the results for the current correlator 
in the full theory, in the p, e and mixed regimes. In 3 we present the results for the partially-quenched 
theory, where the v quarks are quenched, also in the three regimes. Although the e and p regime results 
could have been obtained from earlier literature, we include them for completeness. The mixed-regime 
on the other hand involved a new method to separate the perturbative and non-perturbative modes, that 
will be discussed in detail. In section 4 we present our conclusions and outlook. 

The results in this paper rely heavily on previous similar computations in [8]. We refer the reader to 
those papers for further details of some intermediate steps. 



2. Full Theory Results 

We start by considering SU (Ng + Ny) Yang-Mills theory with Ns flavours with masses nis and A^^ flavours 
with masses rriy. The quark part of the Euclidean continuum Lagrangian reads 

iV„ N 
Le = '^ i^rilixD^ + mv)'^r + X] ^rilixDjj, + ms)lpr , (2.1) 
r=l r=N^+l 

where r is a flavour index; the Dirac matrices are assumed normalised such that 7^ = 7^, {7/^,71/} = 
26 fi^; Djj, is the covariant derivative and N = + Ng. We will consider external sources J that have 
non-zero elements only in the SU^N^) flavour subgroup. 

At large distances, the physics of QCD can be reproduced by chiral perturbation theory. The leading 
order chiral Lagrangian reads 



F2 

-CchPT — ^'^^ 



d^Ud^U^j - -Til^e'^ MU + C/tMe-*^] , (2.2) 



where U G S\J{N). The mass matrix is diagonal with eigenvalues (mt,, m^,, mg, m^) and 9 is the 
vacuum angle. Apart from 6, this Lagrangian contains two parameters, the pseudoscalar decay constant 
F and the chiral condensate S. At NLO in the momentum expansion, additional operators appear in the 
chiral Lagrangian, with the associated low-energy constants Li,L2, ... [13]. For a general N, the number 
of independent couplings is 11-1-2 [14], but for AT = 2 and N = 3, not all of them are independent and 
smaller subsets of 7-1-3 and 10-1-2 couplings respectively are commonly used in this situation [13]. These 
couplings do not depend on the quark masses, but do depend on N. 

One of the simplest correlation functions that can be used to measure F and is also sensitive to the 
NLO couplings Li,L^,LQ and Lg is the left-current two-point function. The numerical advantages of 
such correlator have been discussed in [15]. 
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In QCD, the left-handed flavour current can be formally defined as 

j; = i^T-^^P-^ , (2.3) 

where is a traceless generator of the subgroup SU{Ny), and all colour, flavour, and spinor indices are 
assumed contracted. P_ = (1 — 75)/2 is the left projector. Note that JfJ defined this way is formally 
purely imaginary.^ 

The two-point correlation function between the left-handed currents, averaged over the spatial volume, 
now reads: 

Tr[T-T']C{xo) = J d^x (jS{x)4{0)) . (2.4) 

On the ChPT side, the operator corresponding to Eq. (2.3) becomes, at leading order in the momentum 
expansion, 

= ^IvjT^f/S^i/t] . (2.5) 
The two-point correlation function C{xo) is defined (apart from contact terms) by 

Tr[T''T'']C{xo) = J d'x {jo^{x)J^{0)) . (2.6) 

2.1. The p-regime 

In the p-regime, we express the outcome as a power series in M^/F^, where = 2mE/F'^ is the 
pseudoscalar mass. The power-counting rules for the p-regime are 

M ~ p ~ , (2.7) 

where p is assumed a small quantity, p <C 47r-F. The temporal extent T can in principle be small or large, 
as long as T > 1/ (AttF) . It follows from these assignments that the Goldstone field defined through 
U = exp(2z^/F), behaves effectively as a small quantity, and can be expanded in. Here we have also set 
6 = 0, as is usually done in the p-regime. 

Inserting the Taylor-series of U into Eq. (2.2), the propagator in the quark basis becomes 

U{x) Cdbiy)) = I \Scb6daG{x - y; M^^,) - 6ca6dbE{x - y; M^, , (2.8) 



2 

where 



G(x;M„\)^l^-,^, p^(po,p)^2.(^,^). (2.9) 



*We use this "unphysical" convention since it removes a number of unnecessary overall minus signs from the ChPT 

predictions. 
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V = TL^ is the volume and Af^^ = ^("^'^+"^b) jg mass of a meson constructed out of an a and 6-flavour 
quark, which in practice can be either s or v. On the other hand the singlet contribution is 

ip-x 



with 



+ 



(2.10) 



(2.11) 



[p2 + Ml p2 + Ml 
if Ns + N^^ 0. 

The result for the left-current correlation function in the p-iegime, after spatial integration over the 
source positions and up to contact terms, is: 



{1 + Af) M^P^ixo) - 



N^dG{0;Ml) NsdG{0;Ml) 



F2 dr 



F2 dr 



'"'dM2 



where 



Af = -^g{g-mI)-^g{q-mI) + ^[nsMIu + mI,{ku + l^) 



Am = 



E{Q-Ml,M, 



vv) 



8 

J2 



{N,Ml + N,Ml) (L4 - 2Le) + M2^(L5 - 2Ls) 



The temporal dependence (for \xo\ < T) is contained in the function 

cosh[M^^,(r/2 - \xq\ 



Pv{xo) 



jd?^G{x-Ml) = ^Y. 



^tpoxo 



pI + MI 2M^^sinh[M^„r/2] 



(2.12) 



(2.13) 



(2.14) 



(2.15) 



Up to the second term, which is a constant finite- volume effect, the NLO propagator has the same 
temporal dependence as the LO, if the decay constant and pseudoscalar mass squared are scaled by a 
relative correction given by Ap and Am respectively: 



Fho = F\l + A^), Mlfi^o = + Am) 



(2.16) 



These results agree with those obtained by Gasser and Leutwyler in infinite volume [13] and finite volume 
[1] for Ny = 2 and Ng = 1. Finite volume corrections to F and E have been also obtained to two-loops 
[7]. 

The finite volume corrections can be isolated by 



G(0; M^) = Goo(M^) + Gy(M^) E{0; M\ M'^) = E^{M'\ M'^) + Ev{M'\ M'^) , 



(2.17) 
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where Goo, ^oo are the infinite- volume closed propagators, where instead of momentum sums in eqs. (2.9) 
and (2.10) there are integrals ^. The (finite) functions G\/(M^) and Ev{M'^, M'^) incorporate the volume 
dependence [16]^. 

The limit in which rng ^ ruy, we of course recover the degenerate mass result of [8]. 



2.2. The e-regime 



We consider now the case where all quark masses satisfy m^/^jSF < 1. The results for the correlator 
obtained in a ^-vacuum in [4] and in a fixed-topology in [6, 8] are valid for non-degenerate quark masses: 



2T 



1 + 



where xq = xq/T. The only non-trivial mass dependence is in the function iJ,ai^^'^'"\M): 



U{N) 



YjV r 

PyU + U^Pv (det Uy exp ( — Tr ^MU + M 



(2.18) 



(2.19) 



where Pv is the projector onto the sector of masses rriv, and f^v = niv'^V. 

The constants /3i and koo are related to the (dimensionally regularised) value of 



by 



Introducing p = T/L and 



^(0,0) = 



ap{lo,li) = I dttP-^ 



s{il/t)s'{if/t)- 



(2.20) 



(2.21) 



(2.22) 



where S{x) is an elliptic theta-function, S{x) = X^^_oo 6xp(— vrxn^) = '!?3(0, exp(— ttx)), a numerical 
evaluation of these coefficients is allowed by (see, e.g., Refs. [16, 4]) 



^00 = A ^ 7 XI 



12 4^^sinh2(7rp|n|) ' 
The function /ii(t) appearing in Eq. (3.6) reads (for |r| < 1) 



(2.23) 
(2.24) 



1 

12 



(2.25) 



®The UV divergences of these quantities for d « 4 cancel against those [13] in the Li's as expected. 
®In Ref. [16], the function Gy was denoted by gi. 
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The integral of eq. (2.19) for non-degenerate quark masses can be written in terms of a functional 
derivative: 



m,; 



1 



where 



zf^^^^\Mj)= [ (detC/)'' exp f^TV 
Ju{N) \ 2 



MjU + U^'Mj 



(2.26) 



(2.27) 



where Mj is the block-diagonal matrix: 

Mj = 



{m^ + J)Iv 
m.Js 



(2.28) 



and Iv and Ig are the identity matrices in the v and s subgroups respectively. 
The functional zi^"^'"'' is known in terms of modified Bessel functions [21, 12] 



zi^"^^){Mj) = C,- 



det 



i,j=l,...,N 



(2.29) 



where In{x) is the modified Bessel function and /x| are the eigenvalues of the matrix MjMj multiplied 

by (sy)2. 

In our case we just have to consider two distinct eigenvalues rrig and {my + J)'^. As a useful example 
we consider the case with = 2 and Ng = 1, corresponding to the 2-1-1 flavour QCD: 



+ 



[llyIu+l{llv)Iy{lJ's) - IJ'sIu{lJ'v)Iu+l{l^s)] , 



while 



4^-1) (M). 



(2.30) 



(2.31) 



Another interesting case is that of iVy = 2 and Ng = 0, corresponding to 2 flavour QCD. In this case we 
have 



while 



rf'^^M) 



(2.32) 



(2.33) 
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2.3. The mixed-regime 

Now we turn to the most compUcated case of the mixed regime. In this case, some quarks are in the e 
and some in the p regime and therefore a different factorization of zero and non-zero modes is needed. 
As in the previous sections we start by considering the fuU theory case with the v and s quarks are both 
unquenched. 

The power-counting for this regime is 

rriy ^ e mg p L~ ^ e . (2.34) 

The inspection of the p-regime propagator of eq. (2.8) shows that the modes that become massless in 
the ruy hmit are those corresponding to the generators of SU{Njj). Therefore a factorization that 
would treat the zero-momentum modes of these fields non-perturbatively is 

where Is is the identity matrix in the s sector and Uq G SU{N.u). The perturbative fields ^ satisfy 

d'^x Tr [r'',^] = 0, (2.36) 



I' 



where is a generator of the subgroup SU{Ny). 
It is convenient to include the 6 dependence as 



If the topology is fixed so that 9 is integrated over, the path integral at LO in the e expansion is 

/ df/odet(t7o)^e^^'-[^"(^^°+^»^)], (2.38) 

and the integration over the zero and non-zero modes factorizes. 

The term J(^) is the Jacobian that comes about from the change in the measure when moving from 
[dU] to [dUo] [d^] , which will contribute at NLO as in the e-regime. We describe the computation of this 
measure term in Appendix A. 

The integration over the ^ variables is done in perturbation theory. In order to write the propagator, 
we need to distinguish the indices in the v and s sector, we denote the former by latin letters a, b, ... and 
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the latter by greek ones a, (3,.... The propagator for the ^ fields is : 



^ca{x) Uh{y) 

ica{x)^5b{y) 



Scb^daGix -y;0) - ScaSdb i^E{x - y; 0, 0) - ^^^^^ 



1 



G{x-y;^M, 



ss 



6^pSsaG{x - y; M^) - 5^a5spE{x - y; M^, Mfj 



where G{x, M^) is defined in eq. (2.20) and 



E{x;Ml,Ml)^l,^{l-5l^l) 



"'0;(p2 + M2J(^,2 + M2jF(p)' 



with 



F{p) 



N, 



+ 



(2.39) 
(2.40) 
(2.41) 
(2.42) 

(2.43) 
(2.44) 



.p2 + Ml p2 

The computation of the left-current correlator at NLO, that is at relative order gives a result which 
has the same structure as in the e-regime 



-imixed 



jr2 

{xo) = ^ 

2r2 

+ 



1 



NyG{0,0) + NsG[0 



F2 

/x4°'^'')(M)/ii(xo) 



Ml\ o T2 



SLiNsMl + — (iV^,A:oo + Nsk^o) 



V 



where 



and 



'00 



V 



= T—G 
dT 



/xc7(0'^")(M) = 



2N„ 



(det Uf exp (^Tr 



U + 



(2.45) 



(2.46) 



(2.47) 



lu{N^) 

2.4. Decoupling of the s quarks 

It is useful to rewrite the result of eq. (2.45) in a way which is almost identical to the result in the e-regime 
for a full theory with A^^ degenerate flavours but with a modified F: 



imixed 



\xa) 



El 
2T 

2T 



1 + 



F^ \VV V 

r2 



'00 



2T' 
FW 



l^ai^'^^\M)hi{xo) 



+ ^^00 J 5 



(2.48) 
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where 

1 - ^ (Goo (o, ^ ) - 8L,Ml ) . (2.49) 



P2 _ p2 



The only difference between this expression and that of the full theory with degenerate quarks are 
the finite volume effects in the second line, that are exponentially suppressed in MgsL. 

It is easy to understand these results: the s quarks in the mixed-regime contribute as decoupling 
particles, because the mixed regime probes much lower energy scales than Msg, since the v quarks are 
much lighter and the size of the box is also much larger than the Compton wavelength of the heavy pions: 

Ml^ < < M^^ < (47rF)2. (2.50) 

In this situation one can integrate out the quark within the effective theory [13, 24, 25]. According 
to general symmetry arguments we expect that the theory in this limit can be matched to a theory with 
SU{Ny) flavour symmetry. The effects of the heavy particles can be absorbed in the low-energy couplings 
of the resulting effective theory. Since all the v quarks are degenerate in mass, the result for the correlator 
should be identical to that of eq. (3.6) with = and Ng + Ny ^ Ny which is precisely what we have 
found, apart from exponentially suppressed finite volume effects. In fact the result for the renormalized 
coupling, F^, with Ng = 1 coincides with that obtained in [13] where the matching of the SU{3) flavour 
and the SU(2) flavour effective theories for large strange quark mass was first considered. 

Another observation is that also within the p-regime we can consider a separation of scales < 
Myy <C Mgg. A similar factorization would then be possible for correlators involving only v quarks as 
external legs, and up to exponentially suppressed terms in MggL. The result can be written as the 
correlator in the p-regime for Ny degenerate quarks with mass rriy with modified couplings F as in 
eq. (2.49) and S: 

E ^ E (l - (O, + ^ i|„.„i,,^ . (2.51) 

where 



E'^\x, Ml, Ml) = E^ix, Ml, Ml) - -^Goo(x, Ml), (2.52) 



Ny 

which also coincides with the result of [13]. 

In Figures 1 we show the F and S as functions of M^^/F'^. 

The reason that S does not appear in the mixed-regime of eq. (2.45) is because S appears there only 
at NLO and therefore any correction to it, would be of higher order. 

It is important to stress however that the decoupling only works in a finite volume up to exponentially 
suppressed corrections in MggL, since there is no way to predict these finite volume corrections within 
the effective theory after the s quarks are integrated out. 
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Figure 1: Dependence of the effective couplings F"^ (left) and S (right) in a theory with iVy = 2, and 
with Ng = 1 quarks integrated out, as functions of Mgg/F"^ for F = 90 MeV. The two lines correspond 
to the extreme values of L4,Lq obtained in the phenomenological determinations reviewed in [17]. 

3. Partially- quenched Theory 

We will know consider a partially-quenched theory in which there arc quenched valence quarks of mass 
rriy and Ng sea quarks of mass mg. Note that we consider the generators appearing in the left-currents 
belong to the valence subgroup. 

In order to obtain the partially-quenched results [18], it is simpler to use the replica method of [11]. In 
this method one enlarges the valence sector to Nr degenerate flavours of mass niju , where only A^^ of these 
appear in the external sources. The full symmetry group for zero quark masses is therefore SU{Nr + Ng) 
and the path integral of this theory at the quark level is 

Z[J] = J[dA^] det{p + my + J)^^ det (^ + m^)^"^-^" det(^ + m,)^= g-s^IA^]. (3.1) 

Taking the limit A',. — ^ in this expression, one obtains the supersymmetric formulation of the theory 
[5]. The replica method therefore dictates that one should construct the chiral effective theory for the 
N = Nr + Ng flavours, do perturbative calculations keeping the explicit dependence on and take the 
limit — ^ at the end. We know consider the three regimes in this context. 
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3.1. p-regime 

In the p-regime, as long as Ng 7^ 0, the rephca hmit can be taken and coincides with the result obtained 
using the supersymmetric method [18]. The ^ propagator is that of eq. (2.8) with N^. 
The result for the left-current correlator is that of eq. (2.12) with — > A^^ = and 



E{0; 



G{0;Ml 



^ [NsMlL^ + Mj^Ls) , 



2 M 



F2 



(3.2) 
(3.3) 



In the limit nig rtiy, the full theory result for A^^ degenerate flavours is recovered. Results for the 
meson masses and decay constants at NNLO have been recently obtained [26]. 

It is important to realize that in the partially-quenched theory, the full set of 0{p'^) couplings need to 
be used. The reduction of independent couplings in the full theory with = 2 or = 3 only takes place 
in the unquenched limit, ie. rriy mg. In particular this implies that for Ng = 2, the partially-quenched 
predictions involve more couplings that those that arc physical in the unquenched limit. Obviously these 
couplings cannot be determined from phenomenology (not even in principle) and need to be determined 
on the lattice. 

An interesting observation is that the partially-quenched correction to the meson mass has no 
logarithm in the sea-quark mass, just in the valence quark. If valence quarks masses could be simulated 
in the light regime, for example using Ginsparg- Wilson fermion regularizations [20], and only the sea 
quark masses would be kept unphysically large, the nis dependence of the meson mass would be strictly 
linear at this order of the chiral expansion. In the case of the decay constant the logarithm remains but 
with a smaller coefficient. These features are shown in Figure 2 where we show the dependence of the 
meson mass and decay constant with the sea-quark mass for a value of the valence quark mass of 5 MeV 
for Ng = 2. This is compared with the rUg = rriy dependence in the full theory case, for Ng + = 2. 



In order to recover the fully-quenched case A'^s — ^ 0, it is necessary to keep the singlet meson in the 
theory. When the singlet with a mass iuq is kept in the theory the singlet part of the propagator in 
eq. (2.10) is modified to 



1 



2Nc 



Ei(x- M'^) = — 



(3.4) 



with 



F'^ip) = 1 + 



ap 



2Nr, 



Ng 



+ 



Nr 



p^ + Ml jP + Ml 



(3.5) 



which is well-defined for A^^ = A^^ = 0. Note that as long as either A'^s or A^^ are different from zero, the 
limit mQ — 00 can be safely taken. The results for the two-point function in this limit agree with those 
obtained in [8]. 
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5 10 15 20 5 10 15 20 



Ml Ml 

Figure 2: Chiral correction of (left) and S (right) as a function of mg, with Myy = M^^ = 135MeV in 
the partially-quenched Ng = 2 theory (dashed line) , or of mg = rriy in the full case with A^^. + = 2 (solid 
line), in a lattice of 5 fin. The values of the low-energy couplings have been chosen as the central values 
in the phenomenological determinations of [13, 17]: Li{Mp) = -0.3 • 10"^, Lr,{Mp) = IA-IQ~^ ,L^{Mp) = 
0.9 • 10-3 and Lq{Mp) = -0.2 • IQ-^ with F = 90 MeV. 
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3.2. €-regime 

The results for the partially-quenched theory well all quarks are in the e-regime are 



El 
2T 



1 + 



VVF V 



+ 



2r2 



Pq{N„N^) 



(M)/ii(xo) 



(3.6) 



where 



//a^^(^-^")(M) = lim f 



U{N,+Nr) 



2N„ 



(detU)" exp 



TV 



MU + U^M 



(3.7) 



where is the projector on the SU{Ny) subgroup of SU{Nr). This limit has been defined in [22] as 



-\m) = f 
Jg 



dU J^Tr 



G/(7V,+7V„|iV„) 2A/^, 



PyU + C/^P^J (det exp f^Tr ^MU + U^M 



(3.8) 



where U is an element of the maximal Riemannian manifold, Gl{Ns + A^^|A^^). This integral can be 
obtained as a functional derivative of the functional 



with 



det 



/^i %+i-i(/^i) 



iJ=l,...,Aft, 



nj>i=l,..,Ar„(Mj i^i)W.j>i=N^+l,..,Ni^^j 



{-lYK,{xi) i = l,..,N^ 

Iy{xi) i = N^ + 1,...,2N^ + Ns 



(3.9) 



(3.10) 



lii = jiy.i = I, .■,Ny; Hi = {my + J)T,V, i = N^ + l, ...,2Ny and = //^ = msT^V,i = 2Af^ + l, ...,2Ny + Ns. 

^Pg(iVs,iVt,)^^-j is the same for any value of N^. This is easy to check for small values of and Ng. 
Essentially the choice of N-f, is determined by the dimension of the external source, in eq. (3.9). In 
order to obtain the function iJ,au'^^^"'^'''\M) it is enough to have an external source coupled to one of the 
valence quarks, since due to the SU{Ny) invariance, the contribution of each v quark to eq. (3.8) is the 
same. We can therefore choose the simplest case, A^^ = 1. Any other choice would give the same result. 

Let us consider two simple examples. 

Example 1: = 1, A^s = 1 

It is easy to check that the partition functional is that of a one flavour theory: 



<«(1'1)(M) = /,(m.), 



(3.11) 



while the function iJ,a^'^^^'^\M) is: 



/xa^^(^'i)(M) = /x,S^(Ai.) + ^'f' ii^sL+MlAl^v) - /i/.+i(/x.)/.(M.)] , (3.12) 
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where S'^(/Lt) is the quenched quark condensate in the e- regime at LO[23]: 

S«(/x) = fi [I,{fi)K,{fi) + Iu+i{ti)K,-i{fi)] + -. (3.13) 

A* 

It is easy to check from this expression that the quenched limit is obtained as /x^ — oo: 

hm |Ja^1^^'^\M) = ^i^^1{n^), (3.14) 

and the full theory with just one flavour A^^ = 1 is obtained when the valence and sea masses are the 
same, that is in the limit //s — > /x^: 

hm /xcT^^(^'^)(M) = /xa(°'i)(/x,) = (3.15) 

Note that this last result is a non-perturbative test of Sharpe and Shoresh conjecture that the full theory 
with Ng flavours can be smoothly obtained from the partially-quenched approximation with Ng sea and 
valence quarks. 

Example 2: iV„ = 1, iV^ = 2 

The functional for this case if the one corresponding to a two-flavour theory: 

^P,(2,l)(M) = ^ [I,{^,,f - 7,+i(;,,)/,_i(;x,)] , (3.16) 

while the function /X(T^^*^^'^^(M): 

Ma.^( )(M) = M.S^(^.) - + ^ ^ ^ 



where 

+ —Ii,{fisfKi,+i{fiy)It,+i{fiy) - —I^+i{nsfK^{iJLv)I^{ny). (3.18) 

As in the previous example one can explicitely check that in the limit Ms ^ oo one recovers the quenched 
limit, while in the limit Ms ^ M<^ recovers the full theory with N = 2 degenerate quarks. 

The partially-quenched result then interpolates between the quenched and the full theory with Ns 
flavours. This is shown in left figure of Figure 3 where the function {fiau'^^'^'^\M) — |i/|)//Ut, is shown for 
the partially quenched case as a function of fig and compared with full N = 2 {fig = fi^,) and quenched 
results. On the right figure the /i^, dependence of the condensate for two topological sectors is shown and 
compared with the quenched and full theory, setting //^ = 1. 
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Figure 3: Left: Function {fiGy ' (M) — as a function of Hs = m^ST^ for /Lt„ = niyTiV = 1. The 

star and the dot correspond to the quenched and the fuU (0, 2) resuh respectively at the same yU^,. Right: 
The same function as a function of /n„ for //^ = 1 for topology v = 1 (solid) , compared with the quenched 
function (dotted) and the full (0,2) (dashed). 
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3.3. Mixed-regime 

Both in the p and e regimes we could obtained the partially-quenched result from the full one with 
Ng + Nr quarks by taking the limit N,. ^ at the end of the calculation. While for the current correlator 
the limit seems to be well-defined and indeed is the right answer, in other correlation functions such as 
the pseudoscalar correlator, the limit does not exist. It is easy to see this simply by looking at the ^ 
propagators of eq. (2.39)-(2.42). In the partially-quenched theory the p-regime propagator is the same 
but with Ny — > Nr. If we try to take the iY^ — > limit, the zero-momentum mode contributions in the 
first and third equations explode. This is exactly the same effect that happens in the quenched case if 
the singlet field is not kept in the effective theory [18]. 

In the partially -quenched case, the U{Ns + Nr) singlet can be integrated out and this is true no matter 
whether we are in the p, e or mixed regimes. However, what plays the role here of the non-decoupling 
singlet is the traceless generator of the flavour group SU{Ns + Nr), that is a singlet under the SU{Nr) 
subgroup, whose normalized generator is 



In the mixed- regime not only the SU{Nr) generators become massless in the limit rriy 0, in the limit 
Nr ^ 0, also the pion associated with the rj field gets massless. In fact the propagator for this field can 
be easily derived from eq. (3.27) 



and therefore becomes massless if A^^ ^ and its zero-mode contribution diverges. Note that in the full 
case this is however a massive mode, since the mass goes with the largest massgap. 

The way out of this problem is to modify the factorization in such a way that the zero-momentum 
mode of the rj field is also treated non-perturbatively. So instead of eq. (2.35) we have 



and now the field do not contain the zero-modes of the SU{Nr) generators nor the ij one. As in the 
quenched theory [5], the LO Lagrangian would not factorize into ^ and Uq fields in this case, however it 
does after the integration over 9 is performed, that is in a fixed topological sector. In Appendix B we 
derive the path integral and show that at the LO in the e expansion, the integration over the zero and 
non-zero modes factorize: 




(3.19) 




(3.20) 




(3.21) 




(3.22) 
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The new ^ propagator is: 

[^ca{x)Ub{y) 
'X'ya{x)Csp{y) 



1 
2 

-ScbSsaG [x - y; ^ 



5cbSdaG{x - 2/; 0) - ScaSdbEP^x - y; 0, 0) 
Ml 



1 r 



1 - / Nr n 



SjpSsaGix - y; J 



S^aSsfS ( - y, Ml, Ml) + j^y^ 



where the latin indices refer to the valence and the greek to the sea, and 



E^'^{x-Ml,Ml)^\^Y.{^-^nl 



y {p^ + Ml){p^ + MDFPiip)' 



with 



FPi{p) = 



+ 



Nr 



(3.23) 
(3.24) 

(3.25) 
(3.26) 

(3.27) 
(3.28) 



_p2 + Ml p2 

It is easy to check that the propagator is now well-defined in the limit Nr ^ 0. The singlet part of the 
propagator of the v modes in eq. (3.23) has a double pole structure as in the quenched case, but instead 
of the singlet mass, what appears in the numerator is the heavy mass gap, Ml. This double pole is a 
non-decoupling effect that only appears because the theory is partially-quenched. 

With this parametrization it is easy to check that the left- current propagator is 



+ 



F2 

2T 



1 



F^V^l^' 2 



8L4MI + yk',. 



^a^iiW^M)hi{xo) 



(3.29) 



where fia^'^^^'^"\M) is also the fully quenched result. No double-pole appears for the same reason that 
it did not appear in the quenched case: this observable is not sensitive to them at NLO. 

The decoupling of the s quarks is not possible in the partially-quenched case, because the r] field 
remains light. However, we expect that we should be able to integrate out the scale associated to Mgg 
and match the result to a quenched effective theory. Provided Mgg ^ AttF, this integration can be done 
perturbatively. The quenched Chiral Lagrangian contains additional couplings besides F and S: tjiq and 
a in the standard notation (a mass of the 7] field and a kinetic term). The tree-level matching of m,Q and 
a can be easily read from the propagator of eq. (2.8). The expected p-regime propagator in a quenched 
theory with Ny valence quarks would be of the form 



1 



Uix) idh{y) ) = 7, ScbSdaG{x - y; M^) - 5ca5dbE{x - y; M„^„, M^) 



(3.30) 
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with G = G, and 



E{x;M^^,Ml) = r*- y {^9^ • (3.31) 



Identifying E with £^ in eq. (2.10) with A^^ ^ iV^ = wc find 



a 1 mi Mi, 

° - (3.32) 



2Nc Ns 2Nc Ns 



The quenched F is still the same as that of eq. (2.49), but the matching of the quenched S gives instead 
t (l . ™)) ^ E (l - ^G. (O. f ) . . i|iV.Mi..) . (3.33) 

S gets renormalized in the — limit in the quenched theory, as is well-known. Therefore at one-loop 
the renormalized coupling is 

E, = E(l + mi^)), (3.34) 

while the logs in E^o would cancel on the two sides of the eq. (3.33). The curve S^/S as a function of 
Mgg/F'^ is also shown for A?^s = 1, 2 in Figure 4. 



4. Conclusions 

We have computed the current correlator at next-to-leading order in finite volume Chiral Peturbation 
Theory for non-dcgcncrate quark masses in several interesting regimes. The case when all the quarks are 
in the p-regime could have been easily obtained from earlier literature [13, 18]. However we have also 
considered the case when all quarks are in the e-regime both in the full and partially-quenched theories. In 
the latter case our results provide a non-perturbative check of the Sharpe-Shoresh conjecture, concerning 
the possibility to recover QCD from partially-quenched approximations. 

Finally we have also considered a regime in which some of the quarks masses are in the e and some 
in the p. This mixed-regime required a new zero versus non-zero mode factorization, that we introduced 
and worked out in detail. A further modification was required to treat the partially-quenched theory in 
this case. Our results in the mixed-regime show that the quarks in the p-regime behave essentially as 
decoupling particles, so that the correlator (up to some exponentially suppressed finite volume corrections) 
is that of the e-regime for a theory with a reduced number of flavours (i.e. those in the e-regime), but 
with corrected low-energy couplings by the heavier quarks (i.e. those in the j>- regime). 

These results can be useful for matching lattice QCD and Chiral Perturbation Theory in flnite volumes, 
when the volume is not sufficiently large compared with the Compton wavelength of the lighter pions. 

Clearly the methods developed here can be used for the computation of other correlation functions 
[27]. 
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5 10 15 20 25 

F2 

Figure 4: Dependence of the effective couplings Sr- on M^^/F^ for the partially-quenched theory with 
A^s = 1 (solid) and Ns = 2 (dashed) . F = 90 MeV and the two lines correspond to the two extreme 
values of Lg obtained in the phenomenological determinations reviewed in [17]. Of course in reality -Lg 
can be different in both cases. 
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A. Calculation of the Jacobian 



The parametrization with some or all of the zero modes factorized, that we encounter in the mixed and 

e regimes, has a non trivial Jacobian with respect to the SU{Nf) Haar measure. Here we review first 
how the calculation of the Jacobian can be done for the e regime parametrization at NLO. The method 
is easily extendable to the the mixed regime, so that we briefly mention which are the differences and 
give the results for the mixed regime too. In particular we show that the contribution coming from the 
^ fields f actor izes. 



A metric can be defined through: 

{dUUu) = gabdUadUb 
so that the volume element d/x is obtained as: 

i 

dji = sjdeig jS^ dUi . 



(A.l) 



(A.2) 



The parenthesis (...) means here and in the following that both an integration and a trace are executed 
(/yfTr[...]). 

With the physical fields in the game, our metric matrix g will be of the form 



9 



91 


92 


93 


54 



and we will see that gi contains elements of order 1, g2 = gj contain elements of order e, and 54 
l + 0(e2). 

At next to leading order we are lead by the see saw formula to consider the matrix 
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that has the same eigenvalues and eigenvectors of g to order. 
We will use the parametrization: 



91 








94 - 92 9^^92 



2i4 2ii 

U = UqU^ = e f e f 



(A.3) 



and the expansions: 



im2 



Since ^ only contain nonzero modes, the Z^. satisfy: 

dxJrrn{x) = Q . 

Adding the constant function to the /„. we obtain a complete set: 

XI fmi{x)fmi{y) = 6{x-y)-l . 



(A.4) 



(A.5) 
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The completeness relations for generators such that Tr[r"T*] = ^6"''' read: 



2 

TvlT"- AT"' Fi] = - , 

2N ' ' 2 

1=1 

Ar2-i 

Tr\T"'A]Tr\T"'R] = - , 

2N ' ' ' ' 2 



^ TtIT^AT^B] = -^Tt[AB] + ^Tt[A]Tt[B] (A.6) 

a=l 



Tr[T"'A]Tr[T"'B] = — ^Tr[A]Tr[S] + ^Tr[yl5] (A.7) 

a=l 



and arc valid if A, B are hermitian matrices. 
Inserting (A. 3) in (A.l) we obtain: 



UodUl = - 



2i ; 



1 - e~F^ 



ab 



(A. 



{dUUu) = {uldU^Uodul + dulu^uldU^ + UoduldUoUl + uldulu^dUo) (A.9) 

and note that blocks like U^dU are easy to calculate since they are elements of the algebra. 
The block containing the ^ fields gives: 

uidu^ ^ + - dm + - ^'dc - d^e) (a.io) 

and consequently: 

dUlu^ = {uldU^)^ = -U^dU^ . (A.ll) 
To calculate the block containing the zero modes we define the M{y) matrix by: 

M(y) = e^'^e-^('^+'^<^) (A.12) 

and we see that UodUQ = M{1) — 1. M is a solution for the Cauchy problem: 

^ JA^dne-'^^^'^TT' M(0) = 1 , (A.13) 
dy F 

if F« are the generators of the adjoint irrep ([T",T^] = if'^'T'', [F"']'"' = -if'"'). Finally: 



#"r'' = -A°-^d(t>''T^ = -dUaUl (A.14) 



where (j) = F"'(j)"'. Last equality is a consequence of the reality of the structure constants. 
Exploiting the orthonormality of the /m., for one obtains: 

(gd'^^.lZ = {duyu^uldu^r.,) = A(<5™ + _1_ I ^/„j„,iv[r-2^T-2e - c'r^^r-^]) . (A.15) 

At NLO det(l + a) 1 + Tr[a] if the entries of a arc small. To take the trace of (A.15) one uses the 
completeness relations. The perturbative correction to the determinant is: 

^wo)-i)(-|^)(a • 
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The addend proportional to 6(0) would be there even in the p regime but it is zero in dimensional 
regularization. And this is the explanation why we do not need to consider a measure term in the p 
regime. 

The other correction at NLO comes from Tr[g2 g2] and after straightorward calculations one sees that 
this amounts to 

Tr[glg^%] = j,''^{e). (A.16) 
Combining these results we can calculate the measure: 

dp ~ dp{<P)HaardCJ{0 ^ dp{^) Haard^^ [1 - ^{f}) ■ (A.17) 

The same procedure can be applied to the mixed regime parametrizations. Of course the relations (A. 4) 
and (A. 5) need to be modified properly. We obtained: 



* m€SU{N^) 

^ Jdzjdw J2 Tr[T\CT"-'-T^m^)Tr[T\(T^-T^O]{y) 



b,m€SU{Ny) 

in the factorization for full theory calculations 



JiO c^^^ll-^ [dz Yl TriiT^Y - (r™0'] 

\ meSU{Nv)UTv 

Jdz Jdw Y Tr[T\^T^ - T^i)]{x)Tr[T\(,T^ - T^OMv) 



b,meSU{Ny)UTv 

in the factorization for PQ calculations. 



B. Mixed- regime path integral 



In this appendix we discuss in detail the path integral in the mixed-regime. Starting with the factorization 
in eq. (2.35), the path integral in the absence of sources can be written as 



2. 



d^ J(e) / dUo I de e'"'^ exp (-5) , 

SU{N^) 



(B.l) 



where S can be organized as an expansion in e according to the power-counting of eq. (2.34). The leading 
order is: 



-ps COS 



- -Tr 
2 



+ 



J d^x (rrid^ 



e-rj 



ax)d^ax)] + Micos 



iV., 



TV [p,e] , 



(B.2) 



23 



where the first term if 0{e ^) and the others are all of 0(1). 
Passing to the variables 6 = 6 — 7] and Uq = c^'^^'^Uq: 



J Ju(N^) 

■^^ei^^e-gCos(£)(/..+^(0)^ (B.3) 



where 



^(0 ^ Ml I d'xTr [P^e] + 0{e'), (B.4) 



and B{^, Uq) contains the 0{e^) terms of the original action that do no depend on 6. 
The integral over 6 is therefore of the form: 



/ 



dee^(^\ (B.5) 



where the real part of F(6) has a maximum at ^? = and can be expanded in e via a saddle point 
approximation. Expanding F{6) around this maximum: 



2 V AT J ^ \N, 



F{6) = {]is + A{0) 
the integral can be rewritten as 



+ iiy6, (B.6) 



^ ae^(«) (^1-^ + 0(6^)). (B.7) 
Therefore at the order we need to go, the path integral can be written as: 

/ dUodetiUoreHH^'^^-'^oM)] + B{^, Uo)] . (B.8) 

The LO path-integral factorizes into zero and non-zero mode contributions and the computation of 
correlation functions is then like in the e-regime. 



24 



References 

[1] J. Gasser and H. Leutwyler, Phys. Lett. B 188 (1987) 477; Nucl. Phys. B 307 (1988) 763. 
[2] H. Neuberger, Phys. Rev. Lett. 60 (1988) 889; Nucl. Phys. B 300 (1988) 180. 

[3] P. Hernandez, K. Jansen and L. LcUouch, Phys. Lett. B 469 (1999) 198[arXic:hep-lat/9907022]; 
W. Bietenholz, K. Jansen and S. Shchcredin, JHEP 07 (2003) 033 [hep-lat/0306022]; L. Del Debbio 
and C. Pica, JHEP 02 (2004) 003 [hep-lat/0309145]; L. Giusti, M. Liischer, P. Weisz and H. Wittig, 
JHEP 11 (2003) 023 [hep-lat/0309189]; W. Bietenholz, T. Chiarappa, K. Jansen, K.L Nagai and 
S. Shcheredin, JHEP 02 (2004) 023 [hep-lat/0311012]; L. Giusti, P. Hernandez, M. Laine, P. Weisz 
and H. Wittig, JHEP 01 (2004) 003 [hep-lat/0312012]. L. Giusti, P. Hernandez, M. Laine, P. Weisz 
and H. Wittig, JHEP 04 (2004) 013 [hep-lat/0402002]. T. A. DeGrand and S. Schaefer, Nucl. Phys. 
Proc. Suppl. 140, 296 (2005) [arXiv:hep-lat/0409056]. L. Del Debbio, L. Giusti and C. Pica, Phys. 
Rev. Lett. 94 (2005) 032003 [hep-th/0407052]; H. Pukaya, S. Hashimoto and K. Ogawa, Prog. Theor. 
Phys. 114 (2005) 451 [hep-lat/0504018]; T. Mehen and B.C. Tiburzi, Phys. Rev. D 72 (2005) 014501 
[hep-lat/0505014]; K. Ogawa and S. Hashimoto, Prog. Theor. Phys. 114 (2005) 609 [hcp-lat/0505017]; 
P. H. Damgaard, U. M. Heller, K. SpHttorff and B. Svetitsky Phys. Rev. D 72 (2005) 091501 [hep- 
lat/0508029]. L. Giusti, P. Hernandez, M. Laine, C. Pena, J. Wennekers and H. Wittig, Phys. Rev. 
Lett. 98 (2007) 082003 [arXiv:hep-ph/0607220]. L. Giusti and S. Necco, JHEP 0704 (2007) 090 
[arXiv:hep-lat/0702013] . 

[4] E.G. Hansen, Nucl. Phys. B 345 (1990) 685; E.G. Hansen and H. Leutwyler, Nucl. Phys. B 350 
(1991) 201. 

[5] P.H. Damgaard, M.G. Diamantini, P. Hernandez and K. Jansen, Nucl. Phys. B 629 (2002) 445 
[hep-lat/0112016]. 

[6] P.H. Damgaard, P. Hernandez, K. Jansen, M. Laine and L. Lellouch, Nucl. Phys. B 656 (2003) 226 
[hep-lat/02 11020]. 

[7] G. Golangelo and G. Haefeh, Phys. Lett. B 590 (2004) 258 [hep-lat/0403025]; G. Golangelo, S. Diirr 
and G. Haefeh, Nucl. Phys. B 721 (2005) 136 [hep-lat/0503014]. 

[8] P. Hernandez and M. Laine, JHEP 01 (2003) 063 [hep-lat/0212014]. L. Giusti, P. Hernandez, 
M. Laine, P. Weisz and H. Wittig, JHEP 11 (2004) 016 [hep-lat/0407007]. P. Hernandez and 
M. Laine, JHEP 0610, 069 (2006) [arXiv:hep-lat/0607027]. 

[9] W. Detmold and M.J. Savage, Phys. Lett. B 599 (2004) 32 [hep-lat/0407008]; P.P. Bedaque, 
H.W. Griesshammer and G. Rupak, Phys. Rev. D 71 (2005) 054015 [hep-lat/0407009]. W. Det- 
mold and G.J. Lin, Phys. Rev. D 71 (2005) 054510 [hep-lat/0501007]; 

[10] T. A. DeGrand and S. Schaefer, Phys. Rev. D 72, 054503 (2005) [arXiv:hep-lat/0506021]. A. Hasen- 
fratz, P. Hasenfratz and F. Niedermayer, Phys. Rev. D 72 (2005) 114508 [hep-lat/0506024]. 



25 



H. Fukaya, S. Hashimoto, K. I. Ishikawa, T. Kaneko, H. Matsufuru, T. Onogi and N. Yamada 
[JLQCD Collaboration], Pliys. Rev. D 74 (2006) 094505 [arXiv:liep-lat/0607020]. T. DeGrand, Z. Liu 
and S. Schaefer, Phys. Rev. D 74 (2006) 094504 [Erratum-ibid. D 74 (2006) 099904] [arXiv:hep- 
lat/0608019]. H. Fukaya et al. [JLQCD Collaboration], Phys. Rev. Lett. 98 (2007) 172001 [arXiv:hep- 
lat/0702003]. H. Fukaya et al, arXiv:0705.3322 [hep-lat]. P. Hasenfratz, D. Hicrl, V. Maillart, F. Nie- 
dermayer, A. Schafer, C. Weiermann and M. Weingart, arXiv:0707.0071 [hep-lat]. 

[11] P.H. Damgaard and K. Splittorff, Phys. Rev. D 62 (2000) 054509 [hep-lat/0003017]. 

[12] H. Leutwyler and A. Smilga, Phys. Rev. D 46 (1992) 5607. 

[13] J. Gasser and H. Leutwyler, Annals Phys. 158 (1984) 142; Nucl. Phys. B 250 (1985) 465. 

[14] J. Bijnens, G. Colangelo and G. Ecker, Annals Phys. 280 (2000) 100 [arXiv:hep-ph/9907333]. 

[15] L. Giusti, C. Hoelbling, M. Liischcr and H. Wittig, Comput. Phys. Commun. 153 (2003) 31 [hep- 
lat/0212012]. 

[16] P. Hasenfratz and H. Leutwyler, Nucl. Phys. B 343 (1990) 241. 
[17] G. Ecker, arXiv: hep-ph/0702263. 

[18] C.W. Bernard and M.F.L. Golterman, Phys. Rev. D 46 (1992) 853 [hep-lat/9204007]. S.R. Sharpe, 
Phys. Rev. D 46 (1992) 3146 [hep-lat/9205020]. S. R. Sharpe and N. Shoresh, Phys. Rev. D 62 
(2000) 094503 [arXiv:hep-lat/0006017]; Phys. Rev. D64(2001) 114510 [arXiv:hep-lat/0108003]. 

[19] O. Bar, G. Rupak and N. Shoresh, Phys. Rev. D 67 (2003) 114505 [hep-lat/02 10050]; M. Golterman, 
T. Izubuchi and Y. Shamir, Phys. Rev. D 71 (2005) 114508 [hep-lat/0504013]. 

[20] P.H. Ginsparg and K.G. Wilson, Phys. Rev. D 25 (1982) 2649. D.B. Kaplan, Phys. Lett. B 288 
(1992) 342 [hep-lat/9206013]. Y. Shamir, Nucl. Phys. B 406 (1993) 90 [hep-lat/9303005]; V. Fur- 
man and Y. Shamir, Nucl. Phys. B 439 (1995) 54 [hep-lat/ 9405004]. R. Narayanan and H. Neu- 
berger, Nucl. Phys. B 412 (1994) 574 [hep-lat/9307006]; Nucl. Phys. B 443 (1995) 305 [hep- 
th/9411108]. H. Neuberger, Phys. Lett. B 417 (1998) 141 [hep-lat/9707022]; ibid. 427 (1998) 353 
[hep-lat/9801031]; Phys. Rev. D 57 (1998) 5417 [hep-lat/9710089]. P. Hasenfratz, Nucl. Phys. B 
525 (1998) 401 [hep-lat/9802007] . M. Liischer, Phys. Lett. B 428 (1998) 342 [hep-lat/9802011]. 
Y. Kikukawa and T. Noguchi, hep-lat/9902022. 

[21] R. Brower, P. Rossi and C.I. Tan, Nucl. Phys. B 190 (1981) 699. 

[22] E. Kanzieper, Phys. Rev. Lett. 89 (2002) 250201 [cond-mat/0207745];K. Splittorff and J. J. M. Ver- 
baarschot, Phys. Rev. Lett. 90 (2003) 041601 [arXiv:cond-mat/0209594]. Y. V. Fyodorov and 
G. Akemann, JETP Lett. 77 (2003) 438 [Pisma Zh. Eksp. Teor. Fiz. 77 (2003) 513] [arXivxond- 
mat/0210647]. 



26 



[23] J.C. Osborn, D. Toublan and J.J. Vcrbaarschot, Nucl. Phys. B 540 (1999) 317 [hep-th/9806110]; 
P.H. Damgaard, J.C. Osborn, D. Toublan and J.J. Verbaarschot, Nucl. Phys. B 547 (1999) 305 
[hep-th/9811212]. 

[24] P. Hernandez and M. Laine, JHEP 09 (2004) 018 [hep-ph/0407086]. 

[25] J. Gasser, C. Haefeli, M. A. Ivanov and M. Schmid, arXiv:0706.0955 [hep-ph]. 

[26] J. Bijnens, N. Danielsson and T. A. Lahde, Phys. Rev. D 73 (2006) 074509 [arXiv:hep-lat/0602003]. 
[27] F. Bernardoni, et al in preparation. 



27 



